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The general issue of large N c QCD at nonzero chemical potential is considered with a focus 
on understanding the difference between large N c QCD with an isospin chemical potential and 
large N c QCD with a baryon chemical potential. A simple diagrammatic analysis analogous to 
't Hooft's analysis at \x — implies that the free energy with a given baryon chemical potential is 
equal to the free energy with an isospin chemical potential of the same value plus 1/N C corrections. 
Phenomenologically, these two systems behave quite differently. A scenario to explain this difference 
in light of the diagrammatic analysis is explored. This scenario is based on a phase transition 
associated with pion condensation when the isospin chemical potential exceeds m n /2; associated 
with this transition there is breakdown of the 1/N C expansion — in the pion condensed phase there 
is a distinct 1/N C expansion including a larger set of diagrams. While this scenario is natural, 
there are a number of theoretical issues which at least superficially challenge it. Most of these can 
be accommodated. However, the behavior of quenched QCD which raises a number of apparently 
analogous issues cannot be easily understood completely in terms of an analogous scenario. Thus, 
the overall issue remains open. 
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I. MOTIVATION 

The problem of studying QCD at nonzero chemical po- 
tential is at the crux of nuclear physics. If reliable meth- 
ods could be developed for doing this, then one could 
study cold nuclear matter at zero pressure (which is not 
a bad caricature of the interiors of real nuclei) as well as 
matter under extreme conditions which may be relevant 
in astrophysics or heavy ion reactions. The ultimate pur- 
pose of this paper is to begin to explore the implications 
of the large N c limit of QCD to this problem. Large 
N c QCD has proven to be a useful tool to gain insights 
into hadronic physics and thus it is reasonable to suppose 
that it may prove similarly useful to study QCD with a 
chemical potential 

Before proceeding, it should be noted that the applica- 
bility of large N c QCD to problems in the nuclear physics 
domain such as nuclear matter may be quite problematic 
due to the relevant scales in the problem. In a formal 
large N c expansion one organizes things only by factors 
l/N c . This is clearly useful if N c is large enough and 
if all of the coefficients multiplying the 1/N C factors are 
"natural" . However, if there are a wide variety of scales 
in the problem which have an origin different from 1 /N c 
physics, it may well be that the 1/N C expansion does not 
do a good job of organizing the physics into a hierar- 
chy of scales amenable to expansion. In fact, in nuclear 
physics many important scales are much smaller than 
typical hadronic scales for reasons which have nothing to 
do with large N c physics. For example, the binding of two 
nucleons is formally of order N*. However, the binding 
of the deepest bound two-nucleon state is the deuteron 
at ~ 2 MeV. This is two-order magnitude smaller than 
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the nucleon-Delta mass splitting (~ 300 MeV) which is 
of order 1/N C . Thus one might reasonably worry that the 
vastly different scales in the problem presumably make 
the large N c expansion useless for quantitative or semi- 
quantitative studies of nuclear physics. 

Given this caution, one might conclude that it is point- 
less to use large jV c methods in the nuclear domain. How- 
ever, this would be hasty. In the first place, there may be 
problems for which the scales present in nuclear physics 
do not mix strongly with hadronic scales and which may 
lead to a useful expansion quantitatively. Studies of the 
nucleon-nucleon potential in large N c show that the pat- 
terns of large and small spin-isospin terms in the po- 
tential are consistent with results from straightforward 
1/N C counting 2j. Moreover, even if the large N c expan- 
sion turns out to be of little utility quantitatively it can 
be used to answer questions of principle. For example, 
the question of whether the nucleon-nucleon force is con- 
sistent with a meson-exchange picture can be studied in 
the large N c context yj. The issue is nontrivial since dia- 
grams associated with multi-meson exchange are powers 
of N c larger than the allowed total value. However, can- 
cellations occur between classes of such diagrams to yield 
results consistent with the standard iV c counting. In a 
similar way one might hope to gain qualitative insight 
into QCD at non-zero chemical potential. 

This problem is of practical importance in nuclear 
physics. There is also an underlying theoretical mys- 
tery. At zero temperature the baryon density is strictly 
zero if the absolute value of the chemical potential is less 
than the critical value to make nuclear matter. However, 
from the perspective of the functional integral it is quite 
unclear why the density is zero. This problem has been 
dubbed the Silver Blaze problem^, Q after the Arthur 
Conan Doyle story "Silver Blaze" in which the key clue 
is a dog not barking in the night. Resolving this prob- 
lem is of potential importance since understanding how 
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nuclear matter forms for chemical potential greater than 
this critical value presumably depends on understanding 
how it doesn't form below. It is plausible that large N c 
QCD can give insight into this issue of principle. 

One possible way to gain such insight into the problem 
of QCD at a finite baryon chemical potential is to com- 
pare QCD with an equal chemical potential for up and 
down quarks (which will yield symmetric nuclear mat- 
ter) to QCD with chemical potentials for up and down 
quarks which are equal in magnitude but opposite in sign 
(which makes "isospin matter", i.e., a pion condensate at 
low densities). One reason to focus on this comparison 
is that the isospin matter case is well understood at low 
densities. Firstly at a phcnomenological level this prob- 
lem is amenable to treatment using chiral perturbation 
theory |(J. Secondly, the analogous isospin Silver Blaze 
has been solved: the relevant configurations which con- 
tribute to the functional integral have functional deter- 
minants equal to their zero chemical potential value for 
all chemical potentials less than the critical value (which 
is a chemical potential per quark of m w /2)0- Thus if we 
can understand the difference between the baryon chem- 
ical potential case and the analogous isospin case, we can 
understand the baryon case. At a formal level these two 
systems are quite similar, but at a phenomenological level 
they are very different and our goal is to understand this 
difference. 

On first sight it seems plausible that large N c QCD 
might provide import insights into the differences be- 
tween the two systems. Indeed, a simple diagrammatic 
analysis parallel to the original analysis of 't Hooft (done 
for the zero chemical potential case)0 yields the re- 
sult that the density of up quarks depends only on the 
up quark chemical potential while the density of down 
quarks depends only on the down quark chemical poten- 
tial plus corrections which are suppressed at large N c . 
This implies that baryon matter with a certain chemi- 
cal potential per quark has the same magnitude of quark 
density for up and down quarks as isospin matter of the 
same chemical potential: the only difference being that in 
the isospin case one of the flavors has a negative density. 
A simple way to understand this is that large N c QCD is 
quenched: in quenched QCD quark properties are only 
affected by their propagators and not via the functional 
determinants so there is no way for the chemical potential 
for one flavor of quark to affect the behavior of a quark 
of a different flavor. However, this is at odds with known 
phenomenology. At zero temperature, the critical value 
for the isospin chemical potential goes to zero in the chi- 
ral limit while the critical value of the baryon chemical 
potential presumably does not. Thus, a simple large N c 
analysis of the problem does not resolve the problem of 
interest, rather it yields a new puzzle: namely, why does 
the simple 't Hooft type analysis fail? 

The immediate goal of this paper is to attempt to re- 
solve this problem: before the general behavior of QCD 
at finite baryon chemical potential can be understood, 
this issue must be dealt with. This paper explores a sce- 



nario which appears to provide a natural resolution to 
the problem — the 1/N C expansion breaks down due to 
infrared effects associated with spontaneous breaking of 
U(l) isospin symmetry (associated with the third com- 
ponent) which thereby invalidates the 't Hooft type dia- 
grammatic analysis. Despite its simplicity this explana- 
tion raises some interesting and subtle theoretical issues. 
Indeed, on its face this resolution immediately poses new 
problems. One problem concerns the nature of the break 
down of the 1/N C expansion. A Taylor expansion breaks 
down when it enters a regime in which it ceases to con- 
verge. This does not appear to be the nature of the 
breakdown here — one can be in a regime in which both 
the isospin matter case and the baryon matter case are 
well described by Taylor series in 1/N C and yet they do 
not agree with each other. The fact that they are well de- 
scribed by a converged series suggests that the expansion 
has not broken down, yet the fact that they give distinct 
answers suggests that it has. What is going on? A sec- 
ond problem involves the behavior of the system at small, 
but nonzero temperatures and small chemical potentials. 
In this regime one can use chiral perturbation theory to 
compute the densities and by explicit computation it can 
be seen that the isospin chemical case differs radically 
from the baryon chemical potential case despite the fact 
that the large N c analysis says they should both be equal 
at leading order. This is apparently in contradiction to 
the overall resolution of the paradox: in this regime there 
is no phase transition and hence no apparent reason why 
the 1/N C expansion breaks down. As will be discussed 
in detail in this paper these apparent problems can also 
be resolved and reconciled with the simple explanation 
given above. 

There is an additional problem with the scenario and 
that relates to the quenched approximation to QCD. Lat- 
tice simulations show that quenched QCD incorrectly 
gives the critical point for the baryon chemical poten- 
tial qualitatively (at a chemical potential per quark of 
m 7i-/2)|8j]- Since the diagrammatic analysis implies that 
the large N c limit of QCD is quenched, the known be- 
havior of quenched QCD suggests that whatever problem 
exists at large N c , it is associated with the baryon num- 
ber case and not the isospin case. This seems to be at 
odds with the general resolutions since it is the isospin 
chemical potential which leads to spontaneous symme- 
try breaking. There is a plausible way to evade this 
difficulty — that quenched QCD undergoes a phase tran- 
sition at a chemical potential per quark of associ- 
ated with breaking a symmetry which is not present in 
full QCD. In this picture, quenched QCD with a baryon 
chemical potential behaves analogously to QCD with an 
isospin chemical potential. However, as will be discussed 
in detail this explanation is problematic when viewed at 
a diagrammatic level. 

The difficulties associated with a diagrammatic de- 
scription of quenched QCD raises doubts about the sce- 
nario presented here. It is plausible, however, that the 
problem is associated with peculiarities of quenched QCD 
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(which, after all is not even a physical theory). If this is 
the case, then the explanation outlined here may well be 
correct. 

This paper is organized as follows. In the next section, 
the underlying theoretical differences between the two 
systems (QCD with isospin chemical potential and QCD 
with a baryon chemical potential) is briefly reviewed as is 
the phenomenology of the two cases. In the following sec- 
tion, a simple large N c analysis based on diagrammatics 
is given which yields the prediction that baryon matter 
with a certain chemical potential per quark has the same 
magnitude of quark density for up and down quarks as 
isospin matter of the same chemical potential. This is 
the principal problem addressed here. Next comes a sec- 
tion discussing how a phase transition associated with 
pion condensation appears to resolve this problem and 
various theoretical issues associated with this resolution. 
The final section discusses the challenges to the scenario 
posed by quenched QCD. 

II. BACKGROUND 

For simplicity the problem will be studied in the limit 
of unbroken isospin so the up and down quark masses 
will be taken as degenerate and electromagnetic effects 
ignored. The problem is most easily expressed in terms of 
QCD with a separate chemical potential for up and down 
quarks; that is, consider a system whose Lagrangian is 



given by 

C = £qcd + p u uu + [i d dd . (1) 

With Lagrangian one can evaluate the partition 
function — at least in principle — and from the partition 
function one can compute the free energy density: 

CI,, ,, t\ (log(Z(p u ,p d ,T)) 
y(Pu,Pd,-l) = -j^ (I) 

where V is the spatial volume of the system and (3 = l/T. 
The densities of the two species of quark are given by 

dG &g 

Pu = --h — , Pd = — ■ (3) 

Two important cases can be distinguished — a baryon 
matter (eg. nuclear matter) for which fi u = = p and 
isospin matter (eg., a pion condensate for /i u = —p c i = 
M): 

Z B (p,T) = Z(p,p,T) , Z/(/i,T) = Z(p, -p,T) 

Gb(p-t) = g(p,p,T) , g I (p,T) = g(p,-p,T) (4) 

The focus of this paper is on the difference between 
g B (p,T) and gi(p,T) in the large N c limit of QCD. 

As seen in ref. [9j, both Z B (p,T) and Zj(n,T) can be 
represented as Euclidean space functional integrals: 



Z B (p,T) = J d[A] Y[ det(p+mi) ( det($)+ m - M7o ) ) 2 e~ 5 ™ 

i— heavy 

Z/(/i,T) = J d[A] Y[ det(^)+m l )det( J p+m-^7o) det(p+ m + p-y ) e~ SYM 

i— heavy 

d[A] | dct(^+ m-/j 7 o)| 2 e~ 5 ™ , (5) 

I 



where YM indicates the Yang-Mills action, "heavy" in- 
dicates all quarks heavier than up or down and m is the 
mass of the light quarks. In these expressions the stan- 
dard boundary conditions are imposed: periodic for the 
gluons and anti-periodic for the fermions with an exten- 
sion in the time direction of /3 = 1 /T. Note in the expres- 
sion for Z B the functional determinant is squared (one for 
each flavor) while for the first form for Zi the chemical 
potential for the two flavors is opposite. The second form 
for Zi exploits the identity det(^>+ m + nj )* = det($> 
+ m — (i 70 ) which can easily be proved Q ■ 

Comparing the second form for Z\ with the form for 
Zb in Eq. ©and the fact that det(Ip+mi) is real and pos- 
itive, one sees that the only difference between the two is 
the phase of the functional determinant in the integrand. 



Thus, the entire difference between the phenomenology of 
nuclear matter from isospin matter (a pion condensate) 
stems from this phase. As observed in ref. the fact 
that Re(det($+m- pjo)) 2 < | det(^+ m - ^70 ) | 2 
allows one to deduce in an a priori way that Zi(/i, T) < 
Zb(p, T). The issue of relevance in this paper is whether 
large N c QCD allows us to deduce anything else about 
the relationship between Zi(/j,,T) and Zb((J-,T). 

Before turning to the large N c limit it is useful to re- 
view what is known about the phenomenology of baryon 
matter and isospin matter. Here the focus will be on the 
zero temperature case. There are two reasons for this fo- 
cus. The first is simply because the Silver Blaze problem 
only arises at zero temperature. The second is the anal- 
ysis in this paper will ultimately be done in the large N c 
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limit where one expects the general phenomena of large 
N c continuum reduction to apply [lfi fTTL 111 fl3l jl4|. 

However, as shown in ref. large N c continuum re- 
duction means that at large N c the expectation of an 
operator will be independent of temperature for all tem- 
peratures below some critical value. Thus it is sufficient 
to study the T = case to understand the full low tem- 
perature behavior at large N c . 

First consider the case of isospin matter. At low chem- 
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ical potential one expects that the system remains in 
the vacuum state for chemical potentials less than m^/2, 
which is the critical point since the pion is the state with 
the smallest mass per unit isospin in QCD. At p, — m^/2 
there will be a second-order phase transition to a dilute 
pion condensate whose density then grows with p. The 
system can be analyzed in chiral perturbation theory[(J 
at tree level and the up and down quark densities are 
given by 



t (/x) = -p d (p) = sign (p) 9 (V 2 - ml) 2pft 1 



16^ 4 



1 + 



2 2 

A 2 



(6) 



where A is a typical hadronic scale (for example m p ). 
From Eq. JfJJl it is apparent that if one can neglect the 
higher order terms in the chiral expansion, then there is a 
second-order transition as a function of p with a critical 
p of m,r/2; there is no discontinuity in the density at the 
critical point. On very general grounds the higher-order 
terms are not expected to destroy this qualitative picture: 
there will be a second-order transition at p = m w /2. 

The case of baryon matter is quite different. The tran- 
sition as a function of p is first order at T = 0. Be- 
low a critical value of the chemical potential the den- 
sity is zero. Immediately above this value the den- 
sity is given by the saturation density of nuclear matter 
Po ~ .16 nucleons/fm 3 . The critical chemical potential is 
given by 



Order 


Class of Diagrams 


N 2 C 


• planar graphs with no quark loops. 


Nl 


• planar graphs with one quark loop that 
forms a boundary of the graph. 


N° c 


• graphs with one nonplanar gluon 
and no quark loops. 

• planar graphs with two nonconcentric quark 
loops that form inner boundary of the graph. 



TABLE I: N c counting of leading classes of QCD diagrams 
from 't Hooft counting rules. All diagrams not in these classes 
are lower order in the expansion. The classes are specified by 
topology. Note that there is some freedom in how graphs are 
drawn. A graph is in a class listed in the table if it can be 
drawn in the topology indicated. 



p = (M N - B/A)/N c 



(7) 



where N c , the number of colors, is three for the physical 
world, Mjv is the nucleon mass and B/A is the binding 
energy per nucleon of saturated nuclear matter which is 
estimated to be approximately 16 MeV. These properties 
of baryon matter can be easily inferred from the extrap- 
olation of the properties of finite nuclei [l6|. 

It is apparent that the behavior of the two systems 
is completely different. The order of the transition is 
different as is the critical chemical potential. Perhaps 
more striking than the numerical differences in the crit- 
ical chemical potentials are the qualitative ones: in the 
isospin case the critical chemical potential is of a chiral 
origin and goes to zero as the quark masses go to zero; 
this is not true in the baryon case. Clearly, the single 
sign difference in Eq. leads to profound phenomeno- 
logical differences. In the following sections the issue of 
what large N c QCD tells us about the way a sign plays 
out dynamically will be investigated. 



III. DIAGRAMMATIC ANALYSIS 

The study of large N c QCD began with 't Hooft's sem- 
inal work Q- The key theoretical tool introduced in 
this was a classification of perturbative QCD diagrams 
into classes according to their N c scaling. The assump- 
tions underlying this classification are very simple: the 
quarks are taken to be in the fundamental representa- 
tion of SU(iV c ) while the gluons are in the adjoint; the 

— 1/2 

coupling constant is taken to scale as N c ' while the 
quark mass is taken to scale as iV°. These rules ensure a 
smooth large N c limit for mesonic properties. The anal- 
ysis is greatly simplified by exploiting the fact that color 
content of a color adjoint gluon is essentially that of a 
fundamental-antifundamental bilinear with a color sin- 
glet contribution removed. However, this color singlet 
combination is a 1 /A 2 effect so that at leading order and 
next-to-leading order as far as color counting is concerned 
the gluon acts like a quark-antiquark pair. This is beau- 
tifully encoded via the introduction of 't Hooft's double 
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line representation for the gluons representing the color 
flow (the fundamental and anti-fundamental indices have 
flow in opposite directions). The quarks are represented 
by single lines representing the flow of the color in the 
fundamental representation. Color conservation requires 
that the color entering a vertex must leave it; thus all 
color lines eventually form closed loops. The N c count- 
ing of a diagram is then given by a factor of N c for each 

— 1/2 

closed color loop and a factor of N c ' for every coupling 
constant. 

The results of this analysis are well known and sum- 
marized in Table ITTT1 The issue of what modifications are 
needed to this analysis of the counting rules when a chem- 
ical potential is included is critical here. However, it can 
be dealt with trivially: formally the analysis of the scal- 
ing goes through without modification, provided that the 
chemical potential is of order N®. The reason for this is 
straightforward: the chemical potential enters a diagram 
only via the quark propagator. If the quark mass and the 
chemical potential are each of order N°, then the quark 
propagator itself is also of order N®. Since the inclusion 
of a chemical potential of order unity does not alter the 
N c dependence of the propagator, its inclusion does not 
alter the N c dependence of the diagrams and Table IIIII 
continues to apply. 

In order to use this topological classification of dia- 
grams to deduce the N c scaling properties of physical 
observables, an additional assumption must be made. 
Namely, that the leading N c for some quantity is fixed 
by the N c counting of the lowest-order class of graphs 
which can contribute to the quantity in a perturbative 
expansion. This assumption is nontrivial since a the- 
ory need not be equivalent to the perturbative expansion 
summed to all orders; intrinsically nonperturbative phe- 
nomena can contribute and a priori there is no guarantee 
that nonperturbative contributions will follow the same 
N c scaling laws as the perturbative ones. However, this 
assumption has worked well for the case of zero chemi- 
cal potential in the sense that there are no known cases 
where intrinsically nonperturbative contradictions spoil 
the N c scaling based on diagrammatics. It seems rea- 
sonable as a tentative first step in the study of nonzero 
chemical potentials to make the same assumption. 

The free energy is simply the sum of all connected dia- 
grams (assuming that the resummed perturbative expan- 
sion is equivalent to the theory). The dependence of the 
free energy on the chemical potential requires the con- 
tribution of at least one quark loop. From Table UTTl it 
is apparent that the leading class of diagrams containing 
at least one quark loop are the one-loop diagrams where 
the quark loop bounds the diagram and contributes at 
order N^; two-quark- loop contributions are at least one 
order down in N c . Since the gluons are flavor neutral, 
one-quark-loop diagrams are either loops with up quarks 
or down quarks (the contributions of heavier quark loops 
are neglected since the focus is on the dependence of the 
up and down quark chemical potential.) This implies 
that the /x-dependent part of the free energy density at 



leading order in the expansion is given by 

S(/x«, H d , T) = N c (/( M „, T) + f(jjL d , T)) (1 + 0(1/JV C )) . 

(8) 

The first term on the right-hand side is the sum of the 
one-up-quark loop contributions while the second term 
represents the one-down-quark loops. Note that func- 
tions f(fi) for the up-quark and down-quark contribu- 
tions are identical due to isospin invariance. The overall 
factor of N c follows from the scaling rule seen in Table 
IIIII It is easy to see why Eq. (JHJ) must hold — the large N c 
limit of QCD is quenched. This implies that the contri- 
bution of the two flavors must be independent — the up 
quark density, for example, cannot depend on the down 
quark chemical potential. This independence is encoded 
in Eq. ©. 

The apparent phenomenological consequences of 
Eq. JSJl are quite stark. Combining eqs. JSJ and (0J yields 

&(/*,T) = N c (f(n,T) + f(-n,T))(l + 0(l/N c )) 
g B Q*,T) = 2N c f(fi u ,T)(l + 0(l/N c )) . (9) 

It is easy to see that /(/x,T) must be an even function 
of fx if the 1/N C expansion holds. Qb is given by Eq. 
However, CPT implies that Qb{—^,T) — Qb{^,T) since 
anti-nuclear matter has the same mass as nuclear matter. 
This immediately implies that if 1 /N c corrections can be 
neglected, then f(ji,T) = f(—fi,T). However, if / is 
even, then Eq. © implies that 

g I (fi,T)=g B ( f ,,T)(i + o(i/N c )) . (io) 

Equation (|10f> states that the free energy of baryonic mat- 
ter is identical to that of isospin matter at the same chem- 
ical potential plus 1/N C corrections. It is at the core of 
the paradox: as discussed in Sec. ^ at a phenomenolog- 
ical level the two systems are radically different yet the 
1 /N c expansion suggests that their differences should be 
small. 



IV. A POSSIBLE RESOLUTION 

The resolution to this paradox, presumably either re- 
sides in the theoretical analysis or in the phenomenolog- 
ical description. 

Superficially, the logic underlying the analysis in the 
previous section was identical to that used by 't Hooft 
in the zero chemical potential case where the analysis 
does appear to be reliable. Thus, it is prudent to first 
focus on the possibility that phenomenology as being the 
culprit. A priori this may not seem to be too implausible. 
After all, as discussed in the introduction, the large N c 
limit of QCD can easily fail for nuclear phenomena given 
the very small characteristic scales involved. Thus, one 
might hope that Eq. (|10l) holds formally; this would imply 
that the fundamental problem is simply that the 1/N C 
corrections swamp the leading term. This would make 
nuclear matter very different from isospin matter. This 
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resolution to the paradox amounts to the assertion that 
if we only lived in a truly large N c world (say N c = 1001 
rather than N c = 3), isospin matter and nuclear matter 
at the same chemical potential would have virtually the 
same density. 

However, this phcnomcnological resolution to the para- 
dox is almost certainly not correct. Note that for the 
isospin case there is a second-order transition at /i = 
TOjt/2 for T — 0. For Eq. JTJJ to hold in a large N c world 
the critical chemical potential for the baryon case would 
also have to be ji = and thus tend to zero in the 

chiral limit. For this to occur one would have to make 
nuclear matter which is so deeply bound that the bind- 
ing energy between nucleons cancels the nucleon's mass. 
This is quite far fetched since the nucleon mass appears to 
have nothing to do with the chiral limit. It is even more 
far fetched when one considers the second-order nature of 
the transition. The density at [i just above the transition 
in the isospin case can be made arbitrarily small so that 
Eq. I |l OB implies that the baryon density would also be 
arbitrarily small just above the transition. However, if 
the density in the baryon case is dilute then the nucleon- 
nucleon interactions would be small and thus could not 
cancel the nucleon mass. 

Having disposed of the phenomenological resolution, 
one is forced to conclude that Eq. I|10|l is wrong and, 
hence, so is the theoretical justification underlying it. 
The issue then becomes identifying what is wrong with 
the standard 1/N C analysis of the diagrams, which as we 
have noted is identical to that used by 't Hooft at zero 
chemical potential. In attempting to resolve this the- 
oretically there are a number of issues that appear to 
complicate the analysis. 

The first is that there is no obvious breakdown of 
the 1/N C expansion. If we concentrate on Qb{^,T) and 
Gi([i,T) and imagine writing them both as Taylor series 
in 1 /N c , there is no indication that the series breaks down 
for either quantity. Consider, for example, the regime 
T = 0, ptQ (J, > m„/2, where /j,q is given in Eq. J7J. In 
this regime the C/b(/z,T) is strictly zero. Thus Qb(h,T) 
is accurately described by a Taylor series in 1/N C with 
all terms zero. On the other hand, Qi(n,T) is finite and 
for small enough values of fi and can be arbitrarily 
well described by Eq. (jSJ. Since this regime is supposed 
to be the regime of validity of chiral perturbation the- 
ory, one expects that a systematic series in pion loops is 
possible since the chiral expansion organizes into a loop 
expansion for pion contributions. However, the loop ex- 
pansion for meson is also the l/N c expansion and hence 
one expects that if the chiral expansion is legitimate so is 
the 1/N C expansion. This suggests that, at least in this 
regime, a 1/N C expansion is valid for both quantities of 
interest. This poses a problem since we expect that the 
most plausible theoretical explanation of the difference 
between Qi and Gb is due to a breakdown of the 1/N C 



expansion; if the expansion does not break down, they 
must agree up to l/N c corrections. 

A second issue is associated with the regime of low 
temperatures and small chemical potentials: fi, T <C 
where, at least formally [A, T and are taken to be of 
order N®. In this regime, chiral symmetry implies that 
the system behaves as a gas of weakly interacting pions 
from which one can directly compute the free energy den- 
sity of isospin matter: 

g T Q*,T) = f(T) (1 + 2 cosh ( / i/T))+0(V 2m "/ T ) 

/(t) = T I S iog i 1 - ^ vwt ^ /t ) en) 

Thus we see that the free energy density in this regime, 
though small, is formally of order N®. This should 
be contrasted to the baryon free energy which, in this 
regime, is exponentially small in a l/N c sense Gb{h-, T) ~ 
e~ Nc . Thus, the two free energies continue to differ in 
this regime. The significance of the behavior in this 
regime becomes clear if one studies Gi(n, T) as a function 
of fx at fixed but small T. Equation (|llf> shows explicitly 
that in this regime there is no discontinuity as a function 
of /i all the way down to /x = 0; in the regime under 
consideration there is no phase transition to a pion con- 
densate. However, since the problem of isospin matter 
not matching baryon matter persists in this regime it 
suggests that the phase transition is not the origin of the 
problem. 

These issues appear to suggest that the problem is not 
associated with a breakdown of the 1/N C expansion for 
the free energies, and is not associated with a discontinu- 
ity of the free energy as a function of /x. However, it will 
be argued here that the cause of the problem is with the 
description of the isospin matter system, and is caused 
by a breakdown of the l/N c expansion which, in turn, is 
associated with a phase transition. 

Recall the crux of the problem. A diagrammatic anal- 
ysis at large N c implies that at leading orders the two 
quark flavors do not "talk to" each other: down quark 
observables do not depend on the up quarks and vice 
versa. If they do depend on each other at leading order, 
than the expansion has broken down. This is precisely 
what occurs for the case of isospin matter. To see this, 
let us return to the regime with T — and /i > m n /2 
and fi, rriir -C A is a typical hadronic scale. In this regime 
one can directly compute properties of isospin matter us- 
ing a lowest-order Lagrangian working at tree level. This 
analysis was originally carried out in ref. [|| and was the 
basis for Eq. (jSJ. The important thing here is that in 
this regime there is a pion condensate which can be com- 
puted easily. Expressed in terms of the up and down 
quark fields it is given by: 
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The phase in Eq. 112|l is arbitrary. The presence of an 
arbitrary phase is a signature of spontaneous symmetry 
breaking — in this case the breaking of a U(l) symme- 
try associated with the third component of isospin. Of 
course, the phase is uniquely fixed if by adding an in- 
finitesimal symmetry breaking term to the Lagrangian of 
the system. This symmetry breaking is due to an insta- 
bility in the infrared. This can be seen by looking at the 
chiral susceptibility 



X x = J d x ({imqj 5 T x q)(-x.) {imq~f 5 T x q)(0)) (13) 

where the expression is for Euclidean space. In the bro- 
ken phase fi > , \x IS divergent while it is finite for 
fi < m w . Regardless of the particular value of the phase, 
Eq. (|12fl implies that there is a nonvanishing expecta- 
tion value of an operator which connects up and down 
quarks. This is an unambiguous indication that the up 
quark and down quarks "talk to" each other in isospin 
matter, contrary to the large N c analysis based on dia- 
grams. This shows that the problem does indeed lie with 
the isospin matter and not with baryon matter and is 
associated with the phase transition due to spontaneous 
symmetry breaking. 

It is easy to see at a diagrammatic level why the sym- 
metry breaking leads to a breakdown of the naive 1/N C 
expansion. Consider the \x fr° m Eq. I|13|) which can be 
understood as a quark loop with two insertions of the 
operator imqj^T x q. The fact that \ x diverges in the 
broken phase indicates that even an infinitesimally weak 
perturbation with these quantum numbers has a finite 
effect and that this finite response occurs at leading or- 
der in the 1/N C expansion. Now suppose one were to 
replace the two insertions by a quark-antiquark pair con- 
necting the points as in Fig. IIVI Note that the graph 
with the additional quark-antiquark pair creates an ad- 
ditional quark loop and hence is formally down compared 
to a single quark loop by a factor of N c . However, if the 
quark-antiquark pair carries the quantum numbers of a 
pseudoscalar with the x component of isospin, it acts on 
the original quark loop precisely in the same way as two 
iracpi^Txq insertions and, hence, in the broken phase can 
make a leading order contribution. 

This explanation is both simple and natural. However, 
it appears to be in conflict with the arguments presented 
above. Thus, it is important to understand just how 
these arguments can be evaded. 

The first issue was the lack of an obvious breakdown 




FIG. 1: Feynman graphs with (planar) gluon lines suppressed. 
Left graph is a one-quark-loop graph with two insertions of 
the operator imq^^T^q (represented by filled ellipses). This 
graph corresponds to the chiral susceptibility \x defined in 
Eq. 1131 . Right graph removes these insertions and replaces 
them with a quark-antiquark pair carrying the same quantum 
numbers. 



of the 1/N C expansion for the quantities of interest — 
namely, the free energy densities. In the regime, /io > 
/i > m„/2, both Qj and Gb have well-defined Taylor ex- 
pansions in 1/N C . How then can one argue that a break- 
down in the 1/N C expansion for Qi to cause it to differ 
from Qb as implied by the diagrammatic 1/N C expan- 
sion? 

The answer lies in the possibility that there are multi- 
ple regimes in the problem, each with a valid 1 /N c expan- 
sion but which are not connected with each other due to a 
breakdown of the 1 /N c expansion at the boundary of the 
regimes. Specifically, the scenario is that there is a well- 
defined Taylor expansion in 1/N C for Qi in the regime 
fi < and another well-defined but disconnected ex- 

pansion for fi > TOjr/2. In this scenario, the regime for 
fi < m 7r /2 behaves according to the large N c behavior as 
deduced from the diagrammatic analysis (i.e., Qj = Qb 
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plus possible higher order corrections). In the regime 
where fi > m,/2 (for T — 0), Gi still possesses a 1/N C 
expansion but it is a different expansion than that given 
by the diagrammatic analysis of 't Hooft. This is possible 
only if there is no smooth way to connect the series for 
the two regimes due to a breakdown of the expansion on 
the boundary /! = m^/2. 

As it happens, the series in l/N c for Qi does break 
down at the boundary between these two regimes. Note 
from Eq. JSJl, that the expression contains a smooth func- 
tion of /i times a step function at /it 2 = m 2 /4. Recall, 



however, that the pion mass itself is a function of N c : 
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where the coefficients ; 
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is the nth term in the expan- 
sion. The step function can be formally expanded in N c 
via insertion of the series in Eq. (|14|> and a subsequent 
Taylor expansion. 
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The series clearly breaks down at the point /j 2 



Wo) 
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since the second term is a divergent <5 function. Thus, the 
scenario in which there is a well-defined 1/N C expansion 
in the symmetry broken phase — but a different 1 /N c ex- 
pansion than in the unbroken phase — is viable. 

Next, let us consider the second issue, namely, the be- 
havior of the system in the regime of low temperatures 
and small chemical potentials, i.e. fi, T <C m n . The ap- 
parent problem arises since formally /i, T and m n are of 
order AT?. Equation (|11|) . however, shows that in this 
regime, isospin matter has a finite free energy density as 
N c — > oo while baryon matter has zero free energy den- 
sity. This is in apparent conflict with the large N c expec- 
tation based on diagrammatics seen in Eq. (JSJ. However, 
it is also clear that the problem is not associated with a 
phase transition to a pion condensed phase since no such 
transition occurs in this regime. This, in turn, appears 
to contradict the claim made above that the failure of 
the large N c analysis is due to such a phase transition. 

In fact, there is no contradiction. While it is true that 
isospin matter and baryon matter do differ in this regime, 
the differences — although of relative order 100% — are, in 
fact, sublcading in the 1/N C expansion. The standard 
1/N C counting yields a number density and a free en- 
ergy density of order for chemical potentials of order 
unity. This can be seen straightforwardly from the 't 
Hooft analysis — the leading order is one loop and hence 
proportional to N c . This behavior is seen explicitly, for 
example, in Eq. JHJ): the standard scaling rules / 2 ~ N c , 
m 2 ~ N®, fi ~ N° yield densities and free energy densi- 
ties of order N% when they are non-zero. However, the 
same scaling rules applied to Eq. I|ll[) gives Qj ~ N%. 
Thus, in the regime considered here, both the isospin 
case and the baryon case have free energies which are 
zero at order AT*. They differ from each other at order 
A° (i.e., at the first subleading order) which is exactly 
what one expects from the diagrammatic analysis. 

The behavior seen here in going from zero to finite 



temperature with no effect on the leading order large N c 
behavior illustrates a rather general feature of large N c 
QCD and is related to the general phenomena of large 
N c continuum reduction^ Hi III 111 III El • The un- 
usual feature in the case of the free energy density of 
isospin matter is that at zero temperature the observ- 
able is strictly zero while at non-zero temperature the 
observable is non-zero but still subleading. 

While the explanation given in this section here is 
rather natural there remain some important open ques- 
tions associated with this scenario. Since the essence of 
the original paradox was seen in a diagrammatic anal- 
ysis, a complete resolution of the issue requires an un- 
derstanding of the issue in terms of diagrams. However, 
the resolution to the problem as given above — while very 
plausible — does not make clear precisely which class of 
diagrams play a role at leading order in 1/N C . 

Consider, for example, QCD with an isospin chemi- 
cal potential. The scenario outlined here — that there are 
two distinct 1/N C expansions, one valid for fi < m^j2 and 
one for fi > m^/2 with the expansions breaking down for 
/j = m^/2 — should be describable in terms of diagrams. 
The diagrams which contribute at large N c in the regime 
fi < are presumably the same set which occur for 

fi = 0: one quark loop graphs with planar gluons and 
the quark loop bounding the graph. However, above the 
phase transition our explanation requires a larger set of 
graphs to contribute at leading order. The argument as- 
sociated with Fig. IIVI implies that in the broken phase 
a divergent Xx implies that even infinitesimal perturba- 
tions with the appropriate quantum numbers can lead 
to leading-order response. Additional quark loops can 
couple with these quantum numbers; although the over- 
all strength associated with adding a loop is suppressed 
by 1 /N c they can nevertheless yield leading order results 
due to infrared enhancements. There is, however, an im- 
portant open question associated with this, namely, just 
which class of diagrams does contribute at leading order 
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in this regime? A plausible conjecture is that all graphs 
with planar glue and with an arbitrary number of closed 
quark loops bounding forming an inner boundary of the 
graph (i.e, with no glue inside the graph). This is plausi- 
ble in that it builds in all of the standard 1 /N c rules that 
hold generally except for the quark loop rule (which we 
know to be violated in this regime). 

A parenthetical remark may be useful here: the lead- 
ing order graphs containing quarks for [i = are con- 
ventionally expressed as single quark loops with all of 
the (planar) glue on the inside. However, as Feynman 
graphs these are identical to single quark loops with all 
of the (planar) glue on the outside. The important point 
is simply that they bound the graph. However, if one 
wants to generalize to a situation with multiple quark 
loops, the only topological way to have all of the loops 
bound the graph is to have the glue on the exterior of 
the quark loops. 

It is not immediately apparent, however, if this con- 
jecture is correct. Nor is there an obvious and straight- 
forward way to verify it. Thus, determining whether the 
conjectured class of diagrams for this regime is correct is 
an important open question. 

A second open question concerns the regime of finite 
baryon chemical potential. The explanation given here 
suggests that the standard diagrammatic treatment for 
large N c QCD as done for /i = (planar graphs with a 
single quark loop bounding the graph) remains valid at 
least up to the phase transition where nuclear matter is 
formed (which occurs at /i = fj,o). What class of diagrams 
contributes in the nuclear matter regime? At present we 
do not know. 

Indeed, the general question of how QCD nuclear mat- 
ter behaves at large N c is at present unanswered. We do 
not know, for example, if the transition as a function of 
H at T = is first order (as in the case of N c — 3) or sec- 
ond order. The need to develop reliable theoretical tools 
to study QCD with a baryon chemical potential near the 
transition to nuclear matter is paramount to any attempt 
in understanding nuclear physics from QCD. One may 
hope that large N c QCD will ultimately provide insight 
into this problem, but at present we lack the tools even 
to study large N c QCD in this regime. 

To summarize this section: a natural explanation of 
the different behavior of large N c QCD with a baryon 
chemical potential from large N c QCD with an isospin 
chemical potential is due to a breakdown of the 1/N C 
expansion for QCD with an isospin chemical potential 
at the point where a pion condensation phase transition 
occurs. The two issues raised in this section which sug- 
gested difficulties for such an explanation are apparently 
reconciled with it. However, there remain open questions 
of how the system ought to be described diagrammati- 
cally in the symmetry broken regimes. 



V. IMPLICATIONS FROM QUENCHED QCD 

The behavior of quenched QCD with a baryon chemical 
potential raises issues in many ways quite analogous to 
those of QCD with an isospin chemical potential. Thus, if 
the scenario given in the previous section is correct, one 
would expect a similar scenario to apply for quenched 
QCD. As we will see in this section, the behavior of 
quenched QCD raises challenges to this scenario. 

Recall the essence of the argument of Sec. lIIII that QCD 
is quenched at leading order large N c , at least in a di- 
agrammatic analysis. However, it is long been known 
that quenched QCD at non-zero baryon chemical poten- 
tial gives perverse results. In lattice studies of quenched 
QCD, the critical baryon chemical potential was found 
to be at and not at /io.|5| As \i exceeds there 

is evidence that the chiral condensate changes nonana- 
lytically and that the baryon density discontinuously in- 
creases from zero. Thus a description of baryon matter 
based on quenched QCD with a baryon chemical poten- 
tial appears to be completely wrong: in full QCD the 
chiral condensate remains at its vacuum value and the 
baryon density remains zero until /i = /io- This immedi- 
ately poses a challenge to the explanation in Sec. H VI since 
it suggests that the problem lies with the treatment of 
baryon matter and not with isospin matter. 

To see why note that the qualitative difference be- 
tween quenched QCD and full QCD at non-zero chemi- 
cal potential may be viewed as another symptom of the 
same basic problem seen in the comparison of baryon 
matter and isospin — a diagrammatic large A^ c analysis 
fails. One can consider both full and quenched QCD at 
any N c . Presumably the two descriptions remain qualita- 
tively different at all N c with the transition for quenched 
QCD occurring at fx = and full QCD well above 

this (since there is nothing special about N c — 3 where 
lattice simulations for quenched QCD were done). On 
the other hand, a diagrammatic analysis says that the 
baryon density in large N c QCD should become equiva- 
lent to the density in quenched large iV c QCD since all 
extra quark loops cost a suppression factor of 1/N C . The 
fact that the two differ, seems to be completely analo- 
gous to the fact that baryon matter and isospin matter 
differ despite being equivalent at the diagrammatic limit. 
Clearly, isospin matter plays no role in this problem and 
thus the origin of the difficulty here cannot be a phase 
transition in isospin matter. Since the original problem 
(an isospin chemical potential versus a baryon chemical 
potential) and this one (full versus quenched QCD with a 
baryon chemical potential) have one thing in common — 
full QCD with a baryon chemical potential — it is natural 
to guess that the problem lies there. 

The behavior of quenched QCD might be reconciled 
with our explanation of the difference between isospin 
matter and baryon matter if quenched QCD undergoes 
a phase transition at finite baryon chemical potential in 
a manner analogous to the breaking in isospin matter. 
However, this transition cannot be associated with the 



10 



breaking of a standard baryon number U(l) symmetry, 
since then one would be faced with the questions of why 
the transition occurs at m^/2 and why the transition 
does not occur in full QCD. 

If one views quenched QCD as merely a prescription 
to evaluate functional integrals with the functional deter- 
minant set to unity, it is hard to see how this possibility 
could be realized. However, this prescriptive interpreta- 
tion of quenched QCD has a major drawback: the phys- 
ical quantities of interest are not obtained by differenti- 
ating a generating functional with respect to sources — 
as one does with legitimate quantum field theories. If 
one formulates quenched QCD in terms of a field the- 
ory with a generating functional then one can see how 
such a symmetry breaking could come about. To do this 
one starts with full QCD and then adds a new set of 
particles: for each quark one adds an associated bosonic 
quark (a colored spin-1/2 particle with bosonic commu- 
tation relations) with the same mass and chemical poten- 
tial as the original quarks [T^. These boson quarks, are 
generally denoted as q. Of course, since these fields vio- 
late the spin and statistic theorem, their inclusion leads 
to some clearly unphysical consequences such as states 
with negative norm. Although the theory is physically 
nonsensical in these ways, one can treat it as a prob- 
lem in mathematical physics still computed in the usual 
way. The bosonic quarks can be formally integrated out 
of the functional integral and contribute a factor det(fl 
+ m — /170 ) _1 to the integrand of the gluonic functional 
integral while the corresponding quark contributes a fac- 
tor det(]p+ m — (J, 70 )• The product cancels and thus the 
usual quark fermion determinant is replaced by unity — 
this is standard quenching. 

The key issue here is that the inclusion of bosonic 
quarks implies new symmetries for the system and such 
symmetries can break leading to a phase transition. If 
quenched QCD does undergo a phase transition (associ- 
ated with the breaking of one of these new symmetries) 
when the baryon chemical potential reaches m,r/2, the 
difference between full and quenched QCD at large N c 
appears to be easy to understand. In essence, quenched 
QCD with a baryon chemical potential then behaves 
analogously to full QCD with an isospin chemical po- 
tential: the naive 1/N C expansion breaking down at the 
transition point. It was suggested by Stephanov on the 
basis of a random matrix model study that one expects a 
phase transition associated with the breaking of a sym- 
metry not present in full QCD 01 • Moreover, from first 
principles one expects quenched QCD to undergo a phase 
transition at fj, — m w . The theory has more degrees of 
freedom than the usual QCD — it has boson quarks as 
well as ordinary quarks. Consider an excitation created 
by the current 

J = uijsii + dijsu . (16) 

The lightest state in the spectrum with these quantum 
numbers i.e. the lightest state obtained by J acting on 
the vacuum might be called a "baryo-pion" — it carries 



nonzero baryon number but in many respects acts like 
a pion. In particular, one expects a phase transition at 
the same \x (for a baryon chemical potential) as ordinary 
pions will condense (for an isospin chemical potential). 

To see this, consider the contribution to the suscepti- 
bility for the current J, namely, 

XJ^J d*x (JHx)J(O)) , (17) 

from a one fermion loop diagram in quenched QCD with 
baryon chemical potential fx and with some particular 
arrangement of gluon lines. It is easy to see from the 
Feynman rules that the contribution is identical to the 
contribution to Xxi the chiral susceptibility defined in 
Eq. 1)13(1 at one fermion loop (with an isospin chemical 
potential fj,) and with the same arrangement of gluon 
lines (modulo an overall multiplicative constant). The 
standard 1/N C expansion holds for isospin matter in the 
regime /i < m^/2, so the one fermion loop planar graphs 
should correctly describe Xx f° r l& r g e N c an d (J- > 
This means that in this regime \ J is proportional to x.x- 
Since Xx diverges as ^ — > ra^/2, so does XJ- This type 
of divergence typically signals the onset of a phase tran- 
sition. Physically one might think of such a transition as 
the condensation of "baryo-pions" . 

Thus, quenched QCD does have an instability at a 
(baryon) chemical potential of m^/2 of the sort typically 
associated with a phase transition. This instability is 
presumably associated with the breaking of a symme- 
try which is not present in full QCD. The appearance 
of such a phase transition is what is needed to reconcile 
the difference of behavior of quenched and full QCD with 
baryon chemical potential with our overall explanation. 

While this state of affairs is encouraging, there remains 
a troubling unresolved issue: It provides no way to un- 
derstand the difference between full and quenched QCD 
with a baryon chemical potential at large N c from a di- 
agrammatic perspective. 

Recall that our explanation for the different behavior 
was that quenched QCD undergoes a phase transition at 
fi = m n /2 associated with the breaking of a U(l) symme- 
try which is not present in the full theory. Regardless of 
the phenomenology, we know that the baryon density in 
quenched QCD is generically given diagrammatically as 
the sum of one quark loop graphs with a single insertion 
of the operator q^oq on that loop. Indeed, the essence of 
quenching is precisely that all quark loops not connected 
to external currents do not contribute; they are exactly 
canceled by bosonic quark loops. In the large N c limit, 
only the planar graphs with the quark line bounding the 
graph contribute. 

A potential problem occurs when we consider the 
regime [i < [Iq, where full QCD has not undergone a 
phase transition. Since full QCD is in the same phase 
as at nonzero fj, < fiQ as at /i = 0, one expects that 
the naive 1/N C analysis based on diagrams to hold in 
this regime: at leading order, the baryon density is given 
by planar gluons with a single quark loop which bounds 
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the graph; an insertion of <?7o9 is attached to the quark 
line. Indeed, the scenario outlined in Sec. II VI is based on 
this assumption. However, this is precisely the same set 
of diagrams that occurs in a calculation of the baryon 
density in quenched QCD at large N c regardless of the 
value of /i. Thus, in the regime m,/2 < /x < (j,q, full 
and quenched QCD at large N c appear to have exactly 
the same diagrams contributing, yet to predict different 
densities. The difficulty is that we have no way to under- 
stand the differing behaviors from a diagrammatic point 
of view. 

Note that the situation is qualitatively different from 
our scenario for isospin matter at /i > m w /2. In isospin 
matter the phase transition at fi = implies a new 

1/N C expansion, and associated with this new expansion 
is a new set of diagrams: those with more than a single 
quark loop. However, an analogous behavior cannot hap- 
pen for quenched QCD — by construction, in quenched 
QCD only single quark loop graphs can contribute. In 
this diagrammatic sense, large N c quenched QCD with 
a baryon chemical potential is not analogous to large N c 
QCD with an isospin chemical potential. 

The upshot of this is that the difference between full 
and quenched QCD at large N c and non-zero baryon 
chemical potential is not fully understood. The argument 
analogous to the one used for the isospin case does not 



resolve the issue at a diagrammatic level. Until this issue 
is fully resolved one must view the scenario presented in 
Sec. IIVI as being provisional. While the scenario seems 
quite plausible when viewed in its own right, the lack 
of understanding the apparently analogous problem in 
QCD suggests the possibility of a common explanation 
for both issues which is presently unknown. Of course, 
it is also possible that the scenario in Sec. IIVI is. in fact, 
correct and that there are subtle issues associated with 
quenching which will ultimately explain the puzzle dis- 
cussed in this section. This is, perhaps, not too implau- 
sible. Quenched QCD, after all, is not a true physical 
theory and it is quite conceivable that some peculiar fea- 
ture of it can explain the diagrammatic issue raised in 
this section while leaving the scenario of Sec. IIVI intact. 
Clearly, however, new ideas are needed to resolve this is- 
sue and clarify the nature of large N c QCD with non-zero 
chemical potentials. 
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